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FOREWORD 


This  report  discusses  the  present  status  of  a  program  to  develop  a  procedure  fo 
the  use  of  ultrasonic  techniques  to  nondestructivelv  detect  flaws  in  the  bonding  of 
layered  materials  of  interest  to  the  U.  S.  Navy.  The  program  was  begun  in  1981  by 
the  Materials  Evaluation  Branch  of  the  Materials  Division  at  NSWC  with  the  main 
portion  of  the  program  being  contracted  out  to  Tetra  Tech,  Inc.,  of  Arlington, 
Virginia.  The  theoretical  development  was  done  by  Tetra  Tech  in  Arlington  and  the 
experimental  work  was  done  at  NSWC,  White  Oak  in  collaboration  with  NSWC 
personnel.  Dr.  Ramesh  Shankar  was  the  Tetra  Tech  program  manager  and  principal 
Investigator,  and  was  assisted  by  Dr.  Stephen  S.  Lane  and  Mr.  Thomas  J.  Paradiso. 
Mr.  Clifford  Anderson  of  NSWC  was  the  program  monitor,  and  Mr.  Jeffrey  M.  Warren, 
Ms.  Susan  N.  Vernon,  and  Dr.  P.  L.  Edwards  collaborated  in  the  laboratory 
experiments  and  data  collection. 
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INTRODUCTION 


The  U.  S.  Navy  makes  use  of  multilayered  structures  in  a  number  of  areas;  and 
has  an  interest  in  the  nondestruction  evaluation  of  such  structures.  The  purpose  of 
this  publication  is  to  report  on  the  status  of  a  program  undertaken  to  use  ultra¬ 
sonics  to  nondestructively  detect  flaws  in  the  bonding  of  structures  consisting  of 
alternate  layers  of  bronze  and  rubber  with  epoxy  bonding.  The  purpose  of  the  pro¬ 
gram  was  to  develop  a  generic  model  to  provide  information  as  to  the  best  signal 
processing  techniques  and  instrumentation  to  use  in  the  analysis  of  multilayered 
structures.  This  program  was  undertaken  by  the  Materials  Evaluation  Branch  of  the 
Materials  Division  of  the  Research  and  Technology  Department  in  1981,  and  the  work 
up  to  the  present  time  has  been  done  primarily  through  contracts  with  Tetra  Tech, 
Inc.,  of  Arlington,  Va. 


PROGRAM  PLANS  AND  ACCOMPLISHMENTS 

Tetra  Tech  described1  their  proposed  contribution  to  the  program  in  "PROPOSAL 
TO  DEVELOP  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES"  dated  22  June  1981.  The  report 
states  "The  specific  objectives  of  the  study  are  to; 

•  to  develop  a  computer  simulation  model  to  determine  the 
acoustic  response  at  any  intermediate  layer  of  a  multi¬ 
layered  structure,  for  arbitrary  layer  dimensions,  acoustic 
parameters,  and  driving  source  function, 

•  design  an  experiment,  collect  and  record  ultrasonic  data 
from  rubber-bronze  composite  specimens  containing  flaws  in 
their  intermediate  layers, 

•  process  the  digital  signals  using  spectral  and  cepstral 
methods  to  enhance  flaw  response  and  minimize  reverbera¬ 
tions,  and 

•  utilize  2-D  image  processing  based  software  developed  at 
Tetra  Tech  to  analyze  ultrasonic  C-scan  data  to  enhance  the 
spatial  flaw  response". 


1 Tetra  Tech,  Inc.,  PROPOSAL  TO  DEVELOP  ACOUSTIC  NDE  OF  MULTILAYERED 
COMPOSITES ,  Tetra  Tech  Proposal  P-27191,  22  Jun  1981,  p.  2-3. 
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The  proposal  also  states  that  "The  model  should  also  provide  information  on  the 
effects  of  medium  anisotropies  and  inhomogeneities  on  signal  response." 


The  initial  contract  included  the  following  Tasks,  as  described  in  the  22  June 
1981  proposal^: 


Task  Is  DEVELOPMENT  OF  AN  ACOUSTIC  MODEL  FOR  MULTILAYERED  STRUCTURES 


Subtask  1.1:  Literature  Search 

Subtask  1.2:  Computer  Implementation  of  Acoustic  Models 
Subtask  1.3:  Model  Response  for  Flaws  in  Intermediate  Layers 

Task  2:  DESIGN  OF  EXPERIMENT  FOR  ULTRASONIC  DATA  COLLECTION  FROM 
RUBBER-BRONZE  COMPOSITE  SPECIMENS. 

Subtask  2.1:  Experimental  Verification  of  Model 
Subtask  2.2:  Design  of  Experiment  for  Data  Collection 
Subtask  2.3:  Data  Verification  Preliminary  Analysis 

The  proposal  also  states  that  "based  on  model  results  and  preliminary  data  analysis 
results,  experiments  will  be  designed  to  collect  and  record  ultrasonic  data  from  a 
number  of  "flawed"  and  "unflawed"  specimens.  The  experimental  protocol  will  include 
the  type  of  instrumentation,  driving  force  function,  transducers  to  be  used,  sam¬ 
pling  rate  and  any  prefiltering  that  may  be  required.  .  .  The  data  will  be  collected 
at  NSWC's  facility  in  White  Oak,  Md." 


In  the  work  following  the  initiation  of  the  program,  Tetra  Tech  personnel  have 
developed  the  computer  models,  have  designed  the  experiments  to  test  the  models, 
have  performed  the  experiments  at  NSWC,  White  Oak,  with  the  collaboration  and  assis¬ 
tance  of  NSWC  personnel  and  have  analyzed  the  experimental  results.  A  bronze-epoxy¬ 
rubber  multiple  layer  structure  was  used  to  check  the  model  and  to  determine  suit¬ 
able  signal  processing  methods.  The  signal  processing  methods  included  compensating 
filters,  deconvolution  for  bandwidth  enhancement  and  resolution  improvement,  and  the 
cepstrum  to  separate  the  overlapping  response  from  different  interfaces.  The 
studies  dealt  only  with  pulses  normal  to  the  interfaces,  and  the  bronze  and  rubber 
layer  thickness  were  such  that  the  pulse  travel  time  through  the  bronze  and  rubber 
layers  were  approximately  the  same.  In  the  final  report2 3  for  the  first  year's 
contract  it  was  recommended  that  angled  beams  be  used  for  such  inspection  rather 
than  the  usual  normal  incidence,  pulse-echo  methods. 


The  second  contract  had  as  the  main  heading  on  the  Work  Statement  the 
following:  "Development  and  Analysis  of  Angled  Beam  Methods  for  Multilayered  Struc- 


2Tetra  Tech,  Inc.,  "PROPOSAL",  p.  16-18. 

3Tetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES;  PHASE  1:  ACOUSTIC 
MODEL  DEVELOPMENT  AND  INSPECTION  OF  BRONZE-RUBBER  STRUCTURES,  Tetra  Tech  TC- 
6139, 22  Jul  1982,  p.  65. 
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tures."  The  theoretical  model  developed  under  the  previous  contract  was  to  be 
modified  for  angled  beam  use.  It  was  then  to  be  used  to  predict  optimum  transducer 
placement  for  use  with  rubber-epoxy-bronze  structures.  The  model  predictions  were 
to  be  compared  with  experimentally  determined  values.  The  factors  to  be  compared 
were  mode  conversion  effects,  transmission  and  reflection  coefficients  as  a  function 
of  angle,  attenuation,  multiple  path  interference,  and  phase  changes  and  inter¬ 
ference  due  to  thin  layers.  Experimental  data  from  bronze-rubber-and  other  speci¬ 
men  would  be  used  as  needed  to  check  the  model. 

The  final  report4  for  the  second  contract  period — the  PHASE  2  report —  is  dated 
March  1983.  It  reports  on  the  modification  of  the  theoretical  model  developed  in 
PHASE  1  to  handle  angled  beams.  The  model-simulated  reflection  series  is  compared 
to  experimental  data  collected  using  a  bronze-epoxy-aluminum  structure  for  the 
purpose  of  checking  the  model.  The  angle-beam  model  was  then  used  to  design  a 
procedure  to  inspect  bronze-epoxy-rubber  multilayered  structures  for  "flaw"  and  "no¬ 
flaw"  conditions,  and  the  procedure  was  used  to  collect  data  in  the  laboratory  under 
flaw  and  no-flaw  conditions.  Various  parts  of  this  report  will  be  considered 
below.  The  report  has  a  number  of  errors  in  it  due  to  typographical  mistakes, 
errors  in  the  scales  of  some  figures,  etc.  A  number  of  them  are  listed  with  their 
corrections  in  Appendix  A. 

In  the  first  section  of  the  PHASE  2  report  the  computer  based  model  for  normal 
acoustic  beams  in  multilayered  systems  was  modified  to  be  used  with  off-normal 
beams4.  The  off-normal  beams  are  complicated  by  mode  conversion,  which  leads  to 
their  being  four  reflection  and  four  transmission  coefficients  that  must  be  known 
for  each  interface.  Tetra  Tech  has  taken  this  into  account,  and  has  developed  a 
procedure  to  handle  this  problem.  The  mathematical  procedures  used  are  quite  com¬ 
plicated  and  have  not,  as  of  this  time,  been  completely  verified.  The  model  allows 
these  inputs  for  each  layers  the  layer  thickness,  density,  the  longitudinal  and 
shear  velocities,  the  attenuation,  and  the  mode  (shear  or  longitudinal)  of  propaga¬ 
tion  in  the  layer.  Different  layers  may  have  different  modes  of  propagation.  The 
model  generates  a  complex  reflection  series  as  a  function  of  time. 

In  the  second  section  of  the  PHASE  2  report4,  experimental  data  is  obtained  and 
Interpreted  as  validating  the  following  four  aspects  of  the  model: 

1. )  transit  time  in  each  layer, 

2. )  mode  of  propagation, 

3. )  amplitude  of  interface  response,  and 

4. )  phase  angle  of  interface  response. 

The  multilayered  specimen  used  consisted  of  a  bronze-epoxy-aluminum  specimen 
with  the  bronze  6.35  mm  thick,  the  aluminum  25.40  ram  thick,  and  the  epoxy  bonding 
them  being  nominally  1.00  mm  thick.  Two  5  MHz  transducers  were  used  in  the  pitch- 
catch  mode.  The  crystal  faces  were  about  6.4  mm  by  6.4  mm,  and  they  were  mounted  on 


4Tetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES:  PHASE  2:  ANGLE 
BEAM  VALIDATION  AND  INSPECTION  OF  BRONZE-RUBBER  STRUCTURES.  TETRA  TECH  TC 
6139,  Mar  1983. 
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plastic  wedges  to  produce  a  33°  shear  wave  in  the  bronze  and  a  45°  shear  wave  in  the 
aluminum. 

The  transit  time  through  the  bronze  determined  experimentally  was  6.64  ysec, 
which  agreed  well  with  their  corrected*  theoretical  value  of  6.42  usee.  The  experi¬ 
mental  transit  time  through  the  aluminum  was  determined  as  23.21  usee,  in  good 
agreement  with  their  theoretical  value  of  23.23  usee. 

Although  the  theoretical  and  experimental  values  of  the  transit  time  quoted 
above  are  in  good  agreement,  the  theoretical  values  were  calculated  assuming  the 
bronze  to  be  isotropic,  and  were  based  on  a  shear  velocity  of  2.36  mm/ usee  provided 
to  Tetra  Tech  by  NSWC.  Further  laboratory  studies  at  NSWC  have  shown  that  the  shear 
velocity  in  the  bronze  varies  as  much  as  10%  with  direction,  and  it  may  vary  more 
than  this.  Also,  in  laboratory  studies  at  NSWC  is  was  found  that  at  1  MHz  the 
agreement  between  the  model's  prediction  of  the  transit  times  and  those  found  exper¬ 
imentally  for  bronze  plates  were  not  in  agreement.  These  measurements  are  described 
in  Appendix  B.  The  lack  of  agreement  at  1  MHz  may  be  the  result  of  the  reflected 
beams  being  transposed  relative  to  the  incident  beam  due  to  the  phase  shift  during 
reflection5*6  or  to  the  anisotropy  of  the  bronze,  or  to  a  combination  of  these 
and/or  some  other  effects.  Good  agreement  at  1  MHz  was  found,  however,  between  the 
theoretical  and  experimental  values  for  a  0.261  inch  thick  aluminum  plate,  and  also 
for  a  1.0  inch  thick  brass  plate. 

The  model  developed  by  Tetra  Tech  predicted  a  shear-wave  propagation  mode  for 
the  bronze-aluminum  sample,  and  the  agreement  between  the  model-predicted  transit 
time  and  the  experimental  values  for  the  bronze-aluminum  at  5  MHz  validated  this 
propagation  mode  at  this  frequency,  but  questions  still  remain  in  the  1  MHz  fre¬ 
quency  range. 

The  validation  of  the  signal  magnitude  predicted  by  the  model  was  done  by 
computing  the  power  spectrum  of  the  transmitted  signals  and  then  measuring  the  power 
around  the  center  frequency.  Peak  amplitude  measurements  could  not  be  used  because 
of  phase  shifts  that  occurred  upon  reflection  and  transmission.  The  average  decay 
per  round  trip  in  bronze  was  measured  to  be  -2.9  dB,  whereas  their  theoretical  value 
was  -3.5  dB.  The  equations  used  by  Tetra  Tech  to  calculate  the  reflection  and 
transmission  coefficients  are  given  in  Appendix  C  along  with  another  set  from  a 
different  source.  Tetra  Tech  gave,  via  telephone,  the  calculated  reflection  coef¬ 
ficient  for  a  shear  wave  inQ the  bronze  at  the  bronze-epoxy  interface  at  an  incident 
angle  of  33°  as  0.648eJ6°*5  .  We  were  not  able  to  confirm  this  value  using  either 
set  of  equations  in  Appendix  C. 

Tetra  Tech's  method  used  to  validate  the  signal  phase  given  by  their  model 


*See  Appendix  A 

5Krautkramer ,  J.,  and  Kraut1  ramer,  H.  Ultrasonic  Testing  of  Materials 
(New  York:  Springer-Verlag  ,  1  9)  p.  40. 

6Brekhovskikh ,  L.  M. ,  Waves  in  Layered  Media,  (Academic  Press,  1980),  p. 
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involved  experimentally  determining  the  phase  spectra  of  different  responses  and 
comparing  with  the  model  predictions.  The  first  set  of  equations  in  Appendix  C  were 
used  in  the  model  predictions.  The  responses  considered  were  for  pulses  that  had 
made  n  trips  back  and  forth  across  the  bronze,  and  identified  as  the  (n,0,0) 
response,  where  n  equaled  3,  4,  or  5.  Pulses  were  also  considered  that  had  made 
one,  two  or  three  trips  back  and  forth  through  the  bronze  and  one  trip  back  and 
forth  through  the  epoxy  and  aluminum,  and  they  are  identified  as  the  (1,1,1), 
(2,1,1),  and  (3,1,1)  response,  respectively. 

A  comparison  of  their  theoretical  and  experimental  phase  determinations  are 
given  here.  The  angles  indicated  are  measured  relative  to  the  (3,0,0)  response. 

TABLE  1.  SIGNAL  PRASE  RELATIVE  TO  THE  (3,0,0)  RESPONSE 

Response  Phase  (degrees) 

_ Experiment _ Theory 


(4,0,0) 

-161 

-133 

(5,0,0) 

-322 

-266 

(1,1,1) 

-35 

-64 

(2,1,1) 

+169 

+162 

(3,1,1) 

+8 

+29 

Tetra  Tech's  conclusion  was  that  these  results  verify  their  model;  however,  the 
fact  that  the  phase  angles  reported  for  the  bronze-epoxy  reflection  were  not  consis¬ 
tent  with  the  calculations  made  as  reported  in  Appendix  C  raises  some  questions. 
Also,  the  experimental  measurements  are  discussed  in  Appendix  D,  and  some  questions 
concerning  them  are  raised  there. 

Tetra  Tech  next  considered  the  configuation  that  would  be  optimum  for  locating 
flaws  in  the  lower  epoxy  layer  of  a  bronze-epoxy-rubber-epoxy-bronze  structure.  On 
the  basis  of  their  analysis  the  conclusion  was  reached  that  the  optimum  configur¬ 
ation  for  this  structure  would  involve  a  32°  shear  wave  in  the  upper  bronze  plate 
and  that  the  responses  to  be  observed  were  the  (3, 1,1,1)  response  for  the  no-flaw 
case  and  the  (3,1,1)  response  for  the  flaw  case.  Contributing  to  the  conclusion 
were  the  following  ideas  and  observations: 

a.  The  (3, 1,1,1)  and  (3,1,1)  responses  occur  roughly  midway  between 
the  (3,0,0)  and  (4,0,0)  responses,  and  therefore  interference 
with  other  reflections  is  minimized. 

b.  Both  responses  would  be  large  enough  to  be  detected,  so  that  the 
flaw  or  no-flaw  determination  would  be  made  on  characteristics 
of  a  measured  signal  rather  than  the  absence  of  a  signal,  which 
absence  might  be  due  to  some  other  reason. 


c.  The  phase  difference  between  the  two  signals  was  calculated  to 
be  130°  apart  which  should  allow  for  good  discrimination  between 
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d.  Under  these  conditions  only  a  shear  wave  would  be  propagated 
through  the  bronze,  and  although  both  shear  and  longitudinal 
waves  would  propagate  through  the  rubber  and  epoxy,  the  shear 
wave  would  be  so  highly  attenuated  in  the  rubber  that  it  could 
be  neglected,  and  the  epoxy  layer  being  so  thin  the  shear  waves 
in  it  were  considered  negligible.  Under  these  conditions,  there 
would  be  only  one  important  wave  in  each  layer,  and  therefore, 
the  developed  model  would  be  applicable. 

In  considering  an  operating  frequency  is  was  necessary  to  take  into  account  that 
resolution  is  increased  with  higher  frequencies,  but  so  is  attenuation,  which  can  be 
quite  high  in  rubber.  Taking  these  facts  into  account  an  operating  frequency  of  1 
MHz  was  chosen  for  the  operation. 

Experimental  results  were  obtained  using  two  1  MHz  transducers  on  a  sample 
consisting  of  a  6.35  mm  thick  bronze  plate  bonded  to  a  3.175  mm  thick  rubber 
plate.  A  bond  to  the  bottom  of  the  rubber  was  simulated  by  pressing  a  couplant- 
covered  bronze  plate  against  the  rubber  plate.  The  transducers  were  clamped  in 
place  so  that  their  position  and  coupling  were  unchanged  as  the  lower  bronze  plate 
was  put  into  contact  with  and  removed  from  the  rubber  plate.  The  change  in  the 
waveforms  with  and  without  the  bronze  plate  were  easily  detected  by  the  eye,  which 
means  that  the  signals  vary  depending  on  whether  or  not  there  is  bonding  at  the 
bottom  of  the  rubber  plate,  and  that  an  appropriate  signal  processing  system  may  be 
developed  to  detect  bonding  and  disbonding  at  this  position  in  the  sample.  It 
should  be  noted  however,  that  changing  the  direction  of  the  transducers  resulted  in 
the  pulses  traveling  in  different  directions  in  the  bronze  and  rubber,  and  because 
of  the  anisotropy  of  the  bronze  this  lead  to  a  variation  of  the  received  signal, 
which  can  complicate  the  determination  of  bonding  and  disbonding. 

CONCLUSIONS  AND  RECOMMENDATIONS 

In  the  development  of  an  ultrasonics  procedure  to  detect  and  locate  faults  in 
multi-layered  structures,  the  PHASE  1  undertaking  considered  normal-incident  longi¬ 
tudinal  waves  in  the  pulse-echo  mode,  and  a  computer-based  model  was  developed  to 
predict  the  phases  and  responses  of  reflections  from  such  a  structure.  It  was 
recognized  that  in  some  cases  the  normal- incidence  configuration  would  be  compli¬ 
cated  by  overlapping  pulse  responses,  especially  for  the  case  where  the  transit 
times  through  different  layers  were  approximately  the  same. 

The  PHASE  2  of  the  development  involved  the  modification  of  the  computer  model 
developed  in  PHASE  1  so  that  it  could  be  used  with  off-normal,  or  angled,  ultrasonic 
pulses.  For  the  angled  beam  the  ultrasonic  system  would  operate  in  the  pitch-catch 
mode.  Analysis  showed  that  to  determine  if  a  fault  existed  at  a  given  interface, 
then  the  beam  angle  should  be  selected  so  that  the  transmitting  and  receiving  trans¬ 
ducer  could  be  placed  so  that  the  following  conditions  would  be  met: 

1.  The  reflected  signal  from  the  interface  would  not  be  received  at 
the  same  time  as  other  signals. 

2.  The  signals  from  the  interface  of  Interest  must  be  large  enough 
to  be  detected  both  if  there  is  a  fault  present  or  if  there  is 
not  a  fault  present . 
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3.  The  signal  phase  shift  upon  reflection  where  there  is  a  flaw 
should  differ  from  the  reflection  where  there  is  no  flaw  by  an 
angle  as  close  to  180°  as  possible. 

The  Tetra  Tech  computer  model  was  designed  to  be  able  to  determine  the  incident 
angles  and  transducer  spacings  at  which  the  above  needs  would  best  be  met.  However, 
because  of  the  disagreements  between  the  model  and  experimental  results  as  discussed 
in  Appendix  B,  it  cannot  be  expected  that  the  model  would  give  the  proper  angles  to 
accomplish  the  above  three  objectives.  Tetra  Tech  did  show,  however,  that  the 
difference  between  a  flaw  and  no-flaw  signal — other  factors  being  unchanged — could 
be  detected  from  the  second  epoxy  layer  in  a  bronze-epoxy-rubber-epoxv-bronze 
layered  structure. 

The  next  step  in  the  program  should  be  a  modification  of  the  computer  model  in 
order  to  take  into  account  the  anisotropy  of  materials  and  the  possibility  of  signal 
transposition.  The  model  should  then  be  checked  experimentally  as  reported  in 
Appendix  B,  to  see  if  it  determines  transit  times  correctly.  This  could  turn  out  to 
be  extremely  complicated  since  the  property  directions  in  materials  to  be  tested 
would  not  usually  be  known.  Although  it  appears  rather  unlikely,  procedures  might 
be  developed  in  which  the  existence  of  anisotropy  is  relatively  unimportant  to  the 
end  result.  Also,  if  a  testing  procedure  is  needed  for  a  given  layered  structure 
made  of  specified  materials,  a  testing  procedure  might  be  developed  for  that  partic¬ 
ular  case  taking  into  account  the  specific  properties  of  the  materials  involved. 
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APPENDIX  A 


FINAL  REPORT  TYPING  AND  CORRECTIONS 


Below  are  listed  some  typographical  errors,  errors  in  diagrams,  etc.,  in  the 
final  report^-*  for  the  second  contract,  and  their  corrections. 


Page  2 


R^  and  R£  are  reflection  coefficients,  not 
acoustic  impedances  as  indicated  down  abouc  one- 
third  from  the  top  of  the  page. 


Page  5 


The  equa^ ion 


should  be 


u  »  u  e 
o 


(iu>.x)/v 


Page  10 


In  Section  1.3.3,  second  paragraph,  the  time 
locations  of  the  (1,0,0),  (2,0,0)  and  (3,0,0) 
responses  are  incorrectly  given  as  1.25,  2.50, 
and  3.75  psec,  whereas  they  should  be*  6.42, 
12.84,  and  19.26  usee  using  the  corrected  transit 
time  of  6.42  psec  on  page  16.  The  first  return 
from  the  bottom  of  the  aluminum  occurs  at  29.65 


Page  11 


The  time  scale  of  Figure  1 .4  should  be  increased 
by  a  factor  of  five  in  order  for  the  corrected 
time  locations  (Page  10)  to  be  properly  indicated 
in  the  plot. 


Page  12 


Page  16 


The  time  scale  on  Figure  1.5  should  be  increased 
by  a  factor  of  five. 


The  equation 


'Bronze 


2  x  6.35 
2.2  cos (33) 


6 .58ps 


should  read* 


'Bronze 


2  x  6.35 
2.36  x  cos( 33) 


6.42ps  , 


lTetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES  -  PHASE  2:  ANGLE 
BEAM  VALIDATION  AND  INSPECTION  OF  BRONZE-RUBBER  STRUCTURES,  Tetra  Tech  TC- 
6139,  Mar  1983. 


*  Corrected  values  furnished  by  Tetra  Tech,  Inc. 
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where  Che  shear  velocity  has  been  taken  as  2.36 
mm/ usee. 

Page  19  The  first  sentence  in  the  last  complete  paragraph 

on  this  page  should  read: 

The  average  decay  per  round  trip  is  2.9  db 
which  is  obtained  from  (2.7  +  6.0/2)/2  ■  2.9. 


Page  27 

■ 

The 

table  of 

values : 

V 

_E_ 

(mm/ us) 

vg 

Density 

(gm/cc) 

0 

_E_ 

Thickness 

(mm) 

Bronze 

4.85 

2.20 

8.9 

3.175 

Epoxy 

2.45 

— 

— 

— 

— 

1.000 

Rubber 

2.20 

1.05 

1.37 

10 

6.35 

should 

read* 

V 

(mm/  us) 

V_s_ 

Density 

(gm/cc) 

P 

^s_ 

Thickness 

(mm) 

Bronze 

4.85 

2.36 

8.9 

6.35 

Epoxy 

2.42** 

— 

— 

— 

— 

1.00 

Rubber 

2.20 

1.05 

1.37 

10 

— 

3.175 

Pages  29,31,32  The  time  scales  in  Figures  3.1,  3.2  and  3.3 

should  be  increased  by  a  factor  of  2  so  that 
they  cover  the  range  from  0  to  50  usee. 

The  beam  angles*  for  the  epoxy  and  rubber 
are  also  incorrect  in  the  figure,  as  are 
some  of  the  indicated  velocities.  The  beam 
angles  may  be  calculated  using  Snell's  Law. 


*  Corrected  values  furnished  by  Tetra  Tech. 

**  An  alternate  longitudinal  velocity  for  the  epoxy  of  2.55  mm/^sec  has  also 
been  given. 
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Pages  36, 


37 ,  38 ,  39  The  time  scales*  In  Figures  4.1,  4.2,  4.3, 

and  4.4  should  be  increased  by  a  factor  of  2 
so  that  they  cover  the  range  from  0  to 
20  usee. 

The  beam  angle*  in  the  rubber  should  be  33° 
rather  than  the  indicated  28° . 

Figures  4.1  and  4.2  are  reproduced  on  pages 
vii  and  viii,  and  the  above  corrections 
should  be  made  there  also. 


*  Corrected  values  furnished  by  Tetra  Tech,  Inc 
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APPENDIX  B 

ANGLE  BEAM  MEASUREMENTS  IN  BRONZE  AT  1  MHz 

The  transducer-plate  configuration  is  shown  in  Figure  B-l.  the  transmitter 
beam  passes  through  the  wedge  to  the  interface  with  the  plate  at  the  angle  9,  and 
enters  the  plate  at  the  angle  $.  It  passes  through  the  plate  reflecting  from  the 
top  and  bottom  surfaces.  According  to  the  model  theory  the  angles  9  and  $  satisfy 
Snell's  law,  and  at  the  reflections  within  the  plate  the  angle  of  reflection  is 
equal  to  the  angle  of  incidence  and  the  reflected  wave  begins  where  the  incident 
wave  strikes  the  surface. 

In  the  experimental  arrangement  the  transmitter  is  in  a  fixed  position  and  the 
receiver  is  moved  away  from  it  and  the  ultrasonic  signals  it  receives  are  observed 
on  an  oscilloscope,  and  may  also  be  recorded  for  analysis.  The  distance  between  the 
points  at  which  successive  peaks  are  expected  to  be  received  is  indicated  by  L  in 
Figure  B-l ,  with 

L  -  2D  tan  *  (B-l) 

where  D  is  the  plate  thickness. 

The  transit  time  t  through  the  plate  is  defined  as  the  time  required  for  the 
pulse  to  follow  its  path  from  one  side  of  the  plate,  across  to  the  other  side  of  the 
plate,  and  then  return  to  the  first  side,  so  that 

t  -  L/(V  sin  $)  (B-2) 

where  V  is  the  pulse  velocity. 


Tetra  Tech,  in  its  validation  of  transit  time,  section  2.3.1  of  The  PHASE  2 
report8-1  measured,  using  5  MHz  transducers,  the  time  for  successive  peaking  of  the 
pulses  as  the  receiver  was  moved  away  from  the  transducer  and  obtained  an  average 
transit  time  of  6.64  ysec.  Substituting  Equation  (B-l)  into  Equation  (B-2)  gives 


t 


2D 

V  cos  $ 


(B-3) 


For  the  6.35  mm  thick  bronze  plate,  an  angle  <p  of  33°,  and  a  velocity  of  2.36 
mm/ usee,  the  model  equation  gives  a  transit  time  of  6. 42 ysec. 


B-1Tetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES-  PHASE  2:  ANGLE 
BEAM  VALIDATION  AND  INSPECTION  OF  BRONZE-RUBBER  STRUCTURES,  Tetra  Tech  TC- 
6139,  Mar  1983. 
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Experimental  values  of  transit  times  t  and  distances  L  between  successive 
peaks,  using  1  MHz  transducers,  were  made  by  NSWC  personnel  using  a  10  inch  by  30 
inch  bronze  plate  with  a  thickness  of  6.73  mm.  Some  measurements  were  made  with 
wave  path  along  the  10-inch  direction  (say  the  A  direction)  and  others  with  the  wave 
path  along  the  30  inch  direction  (say  the  B  direction).  The  average  values  obtained 
for  the  various  experimental  runs  are  given  in  Table  B-l . 


Table  B-l.  EXPERIMENTAL  PULSE  MEASUREMENTS  ON  BRONZE  PLATES 


Run 

Direction 

Average 
Transit  Time 

Average 

L 

Velocity 

Angle  o 

1 

A 

9.8  ysec 

24.5  mm 

2.83  mm/ usee 

61.2° 

2 

A 

9.6 

24.0 

2.85 

60.7° 

3 

B 

9.1 

19.4 

2.59 

55.3° 

4 

B 

8.8 

19.2 

2.68 

55.0° 

5 

B 

8.8 

18.5 

2.59 

54.0° 

The  angle  $  given  in  Table  B-l  was  obtained  using  the  experimental  value  of  L  and 
Equation  (B-l).  Using  Snell's  law,  with  the  incident  angle  in  the  wedge  being  39° 
and  the  wave  velocity  of  2.73  mn/usec  in  the  wedge,  the  value  of  for  a  velocity  of 
2.83  mm/usec  calculates  to  be  40.7°  and  for  2.59  mm/ usee  is  36.7°.  The  theoretical 
transit  time  using  Equation  (B-3)  with  a  $  of  33°,  a  velocity  of  2.36  mm/ usee,  and 
the  plate  thickness  of  6.73  mm  is  6.80usec. 

From  the  above  it  is  obvious  that  the  computer  model  does  not  give  the  same 
transit  time  as  was  determined  experimentally.  There  are  other  discrepancies 
also.  The  experimentally  determined  velocities  using  the  angled  beam  model  do  not 
agree  with  the  values  previously  determined,  and  the  experimental  angles  $  do  not 
agree  with  the  values  given  by  Snell's  law. 

It  is  not  known  at  this  time  just  what  causes  this  disagreement.  Further 
velocity  measurements  were  made,  and  it  was  found  that  shear  waves  normal  to  the 
surface  polarized  in  one  direction  had  a  velocity  of  2.35  mm/ usee  and  polarized 
perpendicular  to  that  direction  had  a  velocity  of  2.61  mm/ usee.  There  were  no 
pieces  of  the  bronze  available  for  velocity  measurements  in  other  directions,  and  so 
there  may  be  greater  velocity  variations  than  the  2.35  and  2.61  mm/usec  noted 
above.  Thus,  variations  of  velocity  with  direction  and  polarization  may  be  a  con¬ 
tributing  factor  to  the  above  differences,  as  may  the  transposition  of  the  reflected 
waves  as  discussed  by  Krautkramer  and  Krautkra'mer®-^,  and  by  Brekhovskikh®- 

Because  of  the  disagreements  in  the  velocity  and  angle  measurements  in  the 
bronze,  similar  measurements  were  made  with  a  0.261-inch-thick  aluminum  plate,  and 
the  agreement  there  was  quite  good.  The  velocities  measured  with  the  path  as  in 


^“^Krautkramer ,  j.,  and  Krautkramer,  H.,  Ultrasonic  Testing  of  Materials 
(New  York:  Springer-Verlag  Inc.,  1969)  p.  40. 

®“^Brekhovskikh,  L.  M. ,  Waves  in  Layered  Media,  (Academic  Press,  1980)  p.  104. 
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Figure  B-l  gave  an  average  velocity  of  3.11  mm/ysec  and  an  average  angle  $  of 
48.6°.  Snell's  law  gives  a  value  of  $  equal  to  47.3°,  the  shear  velocities 
measured  for  a  path  normal  to  the  plate  surface  and  with  polarizations  at  0°,  45° 
and  90°  relative  to  the  edge  of  the  plate  were  3.18,  3.19,  and  3.19  mm/ usee.  Thus 
there  appears  to  be  negligible  dependence  of  the  shear  velocity  on  polarization 
direction,  in  contrast  to  the  observations  relative  to  bronze. 
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APPENDIX  C 

REFLECTION  COEFFICIENTS  AT  BRONZE-EPOXY  INTERFACE 

In  a  monthly  report^"^  for  March,  1982,  Tetra  Tech  gave  the  following  informa¬ 
tion  relative  to  the  calculation  of  the  reflection  and  transmission  coefficients  at 
a  bronze-epoxy  interface. 

"In  general,  when  a  plane  wave  is  incident  on  the  boundary 
between  two  solids,  reflected  and  transmitted  waves  of  both  types — 
shear  and  longitudinal — will  be  created.  In  order  to  calculate 
their  amplitudes,  four  equations  are  required.  These  are  provided 
by  the  statements  of  continuity  of  normal  and  transverse  stress,  and 
of  continuity  of  normal  and  transverse  displacement.  Ewing, 

Jardetsky  &  Press^”^  (1957)  have  presented  these  equations,  and 
their  notation  will  be  used  here. 

"The  linearity  of  the  system  makes  it  possible  to  consider 
Incident  longitudinal  waves  separately  from  incident  shear  waves. 

We  wish  to  find  the  following  coefficients: 

R  -  reflection  coefficient  for  P  waves  (longitudinal 
l*  waves)  without  phase  change 

RpS  -  reflection  coefficients  for  P  waves  with  phase 
change  to  S  waves  (shear  waves) 

Tpp  -  transmission  coefficient  for  P  waves  without  phase 
change 

TpS  -  transmission  coefficient  for  P  waves  with  phase 
change  to  S  waves 

The  transmission  and  reflection  coefficients  for  incident  S 

waves  are  similarly  defined. 

Considering  incident  P  waves  first,  the  development  of  Ewing, 

Jardetsky,  and  Press  (EJP)  can  be  used  to  show  that  the  matrix 


C”*Tetra  Tech,  Inc.,  "N60921-81-C-0300.  Monthly  R&D  Status/Performance  Cost 
Report",  15  Mar  1982. 

C"^Ewing,  W.  M. ,  Jarditsky,  W.  S.,  and  Press,  F.,  Elastic  Waves 

in  Layered  Media,  (New  York:  McGraw-Hill  Book  Co.,  1957). 
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equation  relating  the  reflection  and  transmission  efficients  (sic) 
to  the  material  properties  and  to  the  angle  of  incidence  is: 


-1 

0 

n»i3 

m14 

*pp 

m  m 

1 

0 

1 

m23 

m24 

* 

^s 

0 

1 

0 

“33 

m34 

TPP 

1 

0 

-1 

m43 

m44> 

_TPS_ 

0 

(1) 


To  evaluate  the  unknown  matrix  elements  mjj,  let  a,  8,  and  p  be 
respectively,  the  longitudinal  wave  velocity,  transverse  wave  velo¬ 
city,  and  density.  The  angle  of  all  wave  vectors  is  measured  with 
respect  to  the  normal  to  the  interface  and  is  denoted  by  9.  Then 
the  phase  velocity  c  is  the  same  in  all  media  for  all  wave  types, 
and  is  given  by 


c  ■  (wave  velocity)/sin  (9) 
We  also  define  the  quantity  g  such  that 

2. 


8br  *  2^c  pb  ^b  Pb  “  Sr  pr^ 


with  a  similar  definition  for  gfv  obtained  by  interchange  of  sub¬ 
scripts.  The  coefficients  k.  and  are  defined  to  be 


WiF1' 


and 


'.-VT7 


m13 

"  *br  +  pr/pb 

°14 

■  Mbr 

m23 

’  t1  " 

m24 

-  »br 

m_„ 

33 

■  ^  U  -  *brl 

°34 

m13 

°43 

"  *Ar*br 

ra44 

"  m13 

(2) 


(3) 


(4) 


where  i  can  refer  to  either  b  (for  bronze)  or  r  (for  rubber).  With 
this  notation  the  matrix  elements  can  be  written  in  the  following 
form  when  the  wave  is  incident  onto  the  rubber  from  the 
bronze: 


(5) 


C-2 
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Standard  matrix  methods  can  be  used  to  solve  for  the  reflection  and 
transmission  coefficients.  The  matrix  elements  nu..  for  the  case 
where  the  incident  wave  is  in  rubber  can  be  found  from  those  above 
by  interchanging  the  subscripts  r  and  b  wherever  they  appear. 


"Only  a  slight  modification  is  required  to  allow  for  the  case 
where  a  shear  wave  is  incident.  Denoting  the  square  4x4  matrix  on 
the  left  side  of  Equation  (1)  by  M,  we  have  in  this  case 


M 


(6) 


The  elements  of  M  are  the  same  as  before.  They  are  evaluated  in  the 
same  way  as  before,  depending  on  whether  the  Incident  wave  is  in 
bronze  or  rubber." 

In  the  PHASE  2  report  J  the  following  information  and  equations  are  given 
relative  to  calculating  the  reflection  coefficients: 

".  .  .the  potentials  $  and  i|/  are  introduced  whose  derivatives  give 
rise  to  compressional  and  shear  waves  as  follows: 

u  -  3$/3x  -  3Y/3z  1.2 

w  ■  3$/3z  +  3*F/3x 

where  u  and  w  are  displacements  in  x  and  z,  respectively.  If  the 
form 


iux/v 

u  *  u  e 
o 

is  assumed,  the  derivatives  of  equation  1.2  are  equivalent  to  multi¬ 
plication  by  iu/v  . 

The  equations  for  continutiy  of  displacement  and  stress  at  the 
boundary  take  the  form: 


C"3Tetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES  -  PHASE  2: 
ANGLE  BEAM  VALIDATION  AND  INSPECTION  OF  BRONZE-RUBBER  STRICTURES.  Tetra 
Tech  TC-  6139,  Mar  1983. 
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A1  +  \  “  *13^  +  M14  B2 

B1  '  Bi  “  M23  ^  +  M24  B2 

A1  "  \  "  M33  ^  +  M34  B2 

BI  +  3I  =  M43  A2  +  M44  32 


Here  A  refers  to  the  amplitude  of  the  compress ional  potential  and  B 
to  that  of  the  shear  potential.  The  subscript  1  refers  to  the  upper 
layer,  and  2  to  the  lower  layer.  Reflected  waves  (present  in  the 
upper  layer  only)  are  denoted  by  a  prime  sign  (').  The  coefficients 
M.,  .  are  functions  of  the  material  properties  only,  and  are  detailed 
in  [2] . 


"Only  one  incident  wave  at  a  time  need  be  considered  in  the 
upper  layer.  In  the  case  above,  only  a  shear  wave  is  incident 

there.  Then  the  incident  p  wave  amplitude  Aj  is  zero,  and  Bj  is  the 

amplitude  of  the  incident  wave  potential.  We  divide  each  of  the 

Equations  1.3  by  B^,  and  make  the  following  identifications: 


VB1 

»  r 

sp 

VB1 

-  t 

sp 

0 

VB1 

rss 

VB1 

tss 

Then  the  equations  1.3  can  be  written  in  matrix  notation  as 


r  1 

0 

"  -1 

0 

*13 

*U 

■  ^ 

rsp 

1 

0 

1 

*23 

*24 

rss 

0 

1 

0 

*33 

*34 

csp 

1 

0 

• 

-1 

*43 

*44. 

.tss. 

Denoting  the  square  matrix  of  coefficients  on  the  right  by  H,  the 
solution  for  the  desired  reflection  and  transmission  coefficients  is 
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r  « 

rsp 

l-  - 

o 

rss 

-  H-1 

1 

Csp 

0 

css 

m  . 

- 1 

rH 

f 

1.5 


A  complex  matrix  inversion  routine  is  required  to  evaluate  this 
expression  on  a  digital  computer. 

"Finally,  to  get  reflection  coefficients,  the  coefficients 
given  above  must  be  multiplied  by  an  appropriate  ratio  of  veloc¬ 
ities.  Thus,  if  Vs  is  the  shear  velocity  in  layer  1  and  Vp  the 
compressional  velocity  in  layer  2, 


(7) 


and  so  on  in  accordance  with  equation  1.2.” 

According  to  Tetra  Tech,  the  capitalized  transmission  coefficients  in  Equations 
(1)  and  (6)  should  be  lower  case  as  they  are  in  Equations  1.4  and  1.5,  and  the  above 
Equation  (7)  is  then  applicable. 

In  calculating  the  transmission  and  reflection  coefficients  for  the  bronze- 
epoxy  interface,  the  following  parameters  were  used.: 


Material 


Bronze 

Epoxy 


Longitudinal 

Velocity 

4.85  nm/^sec 
2.42 


Shear 

Velocity 

2.36  mm/^sec 
1.05 


Density 

8.89  gm/cm^ 
2.42 


Using  these  values  and  a  33°  incident  angle  in  the  bronze,  the  reflection 
coefficient  for  an  incident  shear  wave  in  the  bronze  at  the  bronze-epoxy  interface 

was  calculated  to  be  0.462e*8®*^  ,  which  does  not  agree  with  the  Tetra  Tech  value 

of  0.648e168*5°. 


It  should  be  noted  that  the  quantities  r  ,  rgs,  etc.  are  reflection  coef¬ 
ficients  for  the  displacement  potentials,  and  the  corresponding  quantities  R  ,  R  , 
etc.,  are  reflection  coefficients  for  the  pressure. 
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Because  of  the  poor  agreement  between  the  above  calculated  values,  the  follow¬ 
ing  equationsc“^  were  used  to  calculate  the  transmission  and  reflection  coefficients 
for  displacement  waves  for  the  case  of  a  shear  wave  incident  in  Medium  1  to  an 
interface  with  Medium  2: 


m 

■  • 

•  •* 

Z  sin 

29t 

-sin  29 

Zf2 

fl 

t 

sin  26 

Z  cos 

29t 

cos  29 

2Z  sin  9t  cos  9 

-2  sin  9  cos 

\t 

r 

cos  29 

-cos 

0 

t 

cos  9 

-Sln  au 

-sin  9 

lr 

T 

-cos  9 

sin 

sin  9 

-cos  9,. 

cos  9. 

R 

sin  9 

t 

it 

It 

- 

.  • 

where 


Also: 


S 

_  r 


T  5ir 

l 


L 


S-j  is  the  amplitude  of  the  incident  shear  wave, 

Sr  is  the  amplitude  of  the  reflected  shear  wave, 

St  is  the  amplitude  of  the  transmitted  shear  wave, 

Lt  is  the  amplitude  of  the  transmitted  longitudinal  wave,  and 

Lr  is  the  amplitude  of  the  reflected  longitudinal  wave. 


V  V 

f  -  [1  +  2(sin20 - —  )] 

1  V 

Pi 


2  2 
V  V  V 

p2  2  '  2  2 

f2  “  U  +  2  (— sin  0  -  —  )] 

2  V  V\ 

1  P2 


PoVo 
2  2 


P1V1  ’ 


C“^Thurston,  R.  N.,  "Wave  Propagation  in  Fluids  and  Normal  Solids,"  Physical 
Acoustics,  Vol.  1,  Part  A,  1964,  p.  83. 
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where  in  medium  1 

Vj  is  the  shear  wave  velocity, 

Vp^  is  the  longitudinal  wave  velocity, 
p^  is  the  medium  1  density, 

0  is  the  shear  wave  incident  and  reflected  angle, 

0£r  is  the  angle  of  the  reflected  longitudinal  wave, 
and  in  medium  2 

V 2  is  the  shear  wave  velocity, 

Vp2  is  the  longitudinal  wave  velocity 
is  the  medium  2  density, 

8  is  the  transmitted  shear  wave  angle,  and 

is  the  transmitted  longitudinal  wave  angle. 

For  the  case  of  the  incident  shear  wave  at  33°  with  the  normal,  as  before,  the 

reflection  coefficient  r  was  calculated  to  be  0.916e“^^*^  ,  which  does  not  agv-v 
with  either  of  the  previously  obtained  values. 
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APPENDIX  D 

EXPERIMENTAL  SIGNAL  PHASE  MEASUREMENTS 

The  signal  phase  measurements  referred  to  in  the  PHASE  2  report*3***  involves  the 
phase  shift  that  the  pulse  signal  undergoes  due  to  the  two  reflections  it  makes  in  a 
trip  back  and  forth  across  the  plate.  In  Figure  D-l,  for  example,  the  signal  phase 
shift  of  the  (4,0,0)  response  relative  to  the  (3,0,0)  response  would  be  the  phase 
shift  due  to  the  reflection  at  the  top  of  the  plate  added  to  that  due  to  the  reflec¬ 
tion  at  the  bottom  of  the  plate;  the  phase  shift  due  to  the  path  distance  is  not 
included. 

Let  the  distance  that  the  pulse  travels  in  getting  to  where  the  (3,0,0) 
response  is  received  be  x^,  and  to  the  (4,0,0)  response  is  received  be  x^.  Consider 
(3,0,0)  and  (4,0,0)  responses  as  shown  in  Figure  D-2  on  a  common  time  scale  t.  In 
order  to  determine  the  (3,0,0)  signal  phase,  the  signal  is  "windowed",  with  t^  as 
the  time  from  the  beginning  of  the  window.  Similar  windowing  is  done  for  the 
(4,0,0)  response,  as  well  as  for  the  other  responses.  The  windowed  signals  were 
about  1  usee  duration,  and  zeros  were  added  to  bring  the  period  up  to  10  usee.  The 
phase  spectra  were  then  computed  using  an  FFT  algorithm. 

In  order  to  get  the  phase  difference  between  two  responses,  the  Tetra  Tech 
procedure  was  to  window  the  responses  by  observing  on  an  oscilloscope  screen  so  that 
the  signals  began  at  the  beginning  of  the  window.  The  Fourier  transforms  were  then 
calculated,  and  the  phase  differences  obtained  from  the  transforms.  To  understand 
under  what  conditions  this  will  give  the  phase  difference  consider  the  following 
analysis. 

The  Fourier  series  for  the  propagating  pulse  may  be  written  as : 

A  « 

F(t)  »  +  l  Cn  cos[na)(t  -  —)  +  *  J  (D-l) 

n-1  n 


where 

a)  is  the  fundamental  angular  frequency,  6.28x10^  rad/sec, 
n  is  the  harmonic  number,  and 

is  the  phase  angle  for  the  nth  harmonic  at  x  equal  zero. 


0-1 Tetra  Tech,  Inc.,  ACOUSTIC  NDE  OF  MULTILAYERED  COMPOSITES  -  PHASE  2: 
ANGLE  BEAM  VALIDATION  OF  BRONZE- RUBBER  STRUCTURES,  Tetra  Tech  TC-6139, 
Mar  1983. 
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Let  the  window  for  the  (3,0,0)  response  begin  at  .  The  signal  at  x3  is  then,  as 
a  function  of  t 


A  *  x. 

F(t)  -  ■=—  +  l  Cn  cos[nw(t  -  -i.)  +  <|>  ] 

n-1  n 


(D-2) 


and  as  a  function  of  t3,  where  t3  -  t  -  t^. 


A  -  x- 

T3(t3)  -  ?(t,  +  t,)  -  ~  4-  I  C^cosJ nw(tA  +  T,  -  ^  ) 


n-1 


3  v 


*3nJ 


CD-35 


Where  equals  <j>n  plus  the  phase  shifts  that  occur  upon  the  surface  reflections. 
The  phase  angle  as  a  function  of  n  at  x  -  x3  is 


K3n* 


(D-4) 


and  this  is  the  phase  angle  at  x3  that  is  given  by  the  FFT  algorithm.  By  a  similar 
procedure,  the  phase  angle  given  by  the  algorithm  at  is 

*4n  ’  ”“<I4  -  T-  >  +  *4n  <D-5) 

n 

The  angle  is  equal  to  plus  the  signal  phase  shift  A4>  due  to  the  reflections 
at  the  surfaces  in  going  from  x3  to  x^,  thus 

x4 

V  “  —  )  +  *3n  +  M  0>-6) 

n 

and 

4»b  ■  *4n  -  *3n  ■  "“<  V  t3  ~  ~  ~ ~  >  +  <°-7> 

n 

Theoretically,  A$  is  independent  of  frequency,  and  for  that  reason  it  has  not  been 
given  a  subsscript  n.  A  basic  assumption  made  by  Tetra  Tech  was  that  by  visual 
observation  they  could  window  the  response  so  that  A^  ■  A$  .  For  this  to  be  the 
case  it  is  necessary  that 


T 


4 


-  0 


(D-8) 


Since  T. ,  t,,  x.,  and  x.  are  constants,  this  means  that  v  must  not  vary  with 
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frequency.  Letting 


and 


AT  -  T,  -  T. 

4 

Ax  -  x.-  x_, 


(D-9) 

(D-10) 


Equation  (D-8)  becomes 


At 

At  -  —  , 
v  * 


( D— 11) 


where  v  is  the  wave  velocity  in  the  medium,  and  At  is  the  time  between  the  beginning 
of  the  two  windows.  If  At  does  not  satisfy  Equation  (D-ll),  and  it  may  be  difficult 
to  begin  the  two  windows  so  that  it  does,  then  Equation  (D-8)  is  not  satisfied,  and 
setting 


X4  "  *3 

T.  -  T, - - - -  -  C 

4  3  v 


(D-12) 


then  equation  (D-7)  becomes 


Atn  "  Cnw  +  A<j>, 


(D-13) 


in  which  case  the  measured  signal  phase  shift  is  a  linear  function  of  frequency. 

Tetra  Tech  has  assumed  that  they  could  select  t,  and  t,  such  that  C  »  0,  but 
their  data  from  the  PHASE  2  report®"*  which  is  shown  nere  appears  to  have  a  linear 
variation  with  frequency. In  all  the  three  cases  of  TABLES  D-l  and  D-2,  the  phase 
difference  is  due  to  one  transit  path  across  and  back  through  the  bronze. 

Using  linear  regression  to  determine  the  curves,  the  linear  equations  for  the 
three  cases  become: 

for  (3,0,0)  and  (4,0,0)  responses, 

A*  -  8.5  f  -  190.8°  (D-14) 

n 

for  (4,0,0)  and  (5,0,0)  responses, 

A*  -  12  f  -  230°  (D-15) 

n 
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for  (1,1,1)  and  (2,1,1) 

A*  -  18  f  -  241.2°,  (D-16) 

n 

where  is  measured  in  degrees  and  f  is  in  megahertz.  The  coefficient  of  cor¬ 
relation  is  0.986  or  better  for  each  of  the  equations.  The  units  of  the  coef¬ 
ficients  of  f  are  degrees  per  megahertz,  and  to  convert  the  units  to  seconds  as  is 
needed  for  use  in  Equation  (D-13)  they  must  be  multiplied  by  2.78  XI Q“*  sec  MHz  per 
degree. 


TABLE  D-l.  PHASE  DIFFERENCES  FOR  (3,0,0)  -  (4,0,0)  AND  (4,0,0)  - 

(5,0,0)  RESPONSES 

Phase  Difference  (Degrees) 


(3,0,0) 

(4,0,0) 

and 

and 

Frequency  (MHz) 

(4,0,0) 

(5,0,0) 

4.8 

-150 

-172 

5.0 

-148 

-170 

5.2 

-147 

-168 

5.4 

-145 

-166 

5.6 

-143 

-162 

Average  Difference: 

-147±3 

-168±4 

TABLE  D-2. 

PHASE  DIFFERENCES  FOR  (1,1 

,1)  -  (2,1,1)  RESPONSES 

Phase 

Difference 

Frequency 

(1,1.1) 

and  (2,1,1) 

3.6 

-177 

3.8 

-172 

4.0 

-169 

4.2 

-166 

4.4 

-162 

Average  Difference 

-169  ± 
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The  above  equations  then  become 


A*  -  2.36  X  10“8nw  -  190.8°  (D-17) 

n 

Ai|i  -  3.33  x  10  "8nu>  -  230°  (D-18) 

n 

Ai|)  -  5.00  x  10“8nai  -  241.2°  (D-19) 

n 

A  comparison  of  these  equations  with  Equation  (D-13)  shows  that  the  constant 
term  of  the  equations  is  just  A$  of  Equation  (D-13),  which  is  just  the  signal  phase 
shift  in  going  from  the  first  pulse  under  consideration  to  the  second  one.  They  may 
be  compared  with  Tetra  Tech's  experimental  and  theoretical  values  shown  in  Table 
1.  There  is  disagreement  in  both  cases.  The  disagreement  with  the  experimental 
value  results  from  the  assumption  of  an  error  in  windowing,  and  the  disagreement 
with  the  theoretical  value  may  be  due  to  the  reflection  calculations  considered  in 
Appendix  C. 

Comparison  of  the  equations  with  equation  (D-13)  also  shows  that  C  is  equal  to 
the  multiplier  of  nu>  in  the  equations.  From  Equations  (D-8)  and  (D-12)  it  is  seen 
that  a  non-zero  C  is  due  to,  and  equal  to,  errors  in  the  values  of  At,  which  results 
from  errors  in  windowing;  that  is,  in  errors  in  the  beginning  times  of  the 
windows.  From  the  Equations  (D-17),  (D-18),  and  (D-19)  it  is  seen  that  windowing 
errors  of  the  responses  considered  vary  from  0.0236  to  0.0500  microseconds.  This 
implies  an  error  in  the  timing  of  the  beginning  of  the  window  that  is  2Z  to  5%  of 
the  window  duration.  In  looking  at  the  windowed  pulses  on  pages  17  and  18  of  the 
PHASE  2  report,  it  appears  that  such  errors  could  easily  be  made. 
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